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Abstract 

Epitaxial self-assembled quantum dots (SAQDs) represent an important step in the advancement of semi- 
conductor fabrication at the nanoscale that will allow breakthroughs in electronics and optoelectronics. In 
these applications, order is a key factor. Here, the role of crystal anisotropy in promoting order during 
early stages of SAQD formation is studied through a linear analysis of a commonly used surface evolution 
model. Elastic anisotropy is used a specific example. It is found that there are two relevant and predictable 
correlation lengths. One of them is related to crystal anisotropy and is crucial for determining SAQD order 
Furthermore, if a wetting potential is included in the model, it is found that SAQD order is enhanced when 
the deposited film is allowed to evolve at heights near the critical surface height for three-dimensional film 
growth. 



1 



Epitaxial self-assembled quantum dots (SAQDs) represent an important step in the advance- 
ment of semiconductor fabrication at the nanoscale that will allow breakthroughs in electronics 
and optoelectronics. SAQDs are the result of linearly unstable film growth in strained heteroepi- 
taxial systems such as SixGei_a;/Si and In^.Gai_^As/GaAs and other systems. SAQDs have great 
potential for electronic and optoelectronic applications. In these applications, order is a key factor. 
There are two types of order, spatial and size. Spatial order refers to the regularity of SAQD dot 
placement, and it is necessary for nano-circuitry applications. Size order refers to the uniformity 
of SAQD size which determines the voltage and/or energy level quantization of SAQDs. It has 
been observed that crystal anisotropy can have a beneficial effect on SAQD order. ' Here, the role 
of crystal anisotropy in promoting order during early stages of SAQD formation is studied through 
a linear analysis of a commonly used surface evolution model . ^'^'^'"^^^ 

The linear analysis addresses the initial stages of SAQD formation when the nominally flat film 
surface becomes unstable and transitions to three-dimensional growth. This early stage of SAQD 
growth determines the initial seeds of order or disorder in an SAQD array, and can be analyzed 
analytically. A dispersion relation as in^'^ is used, and two predictable correlation lengths that 
grow as the square-root of time (Eqs.[l4]and[l7]) are found. One of the correlation lengths results 
from crystal anisotropy. This length plays a limiting role in the initial order of SAQD arrays, and it 
is shown that anisotropy is crucial for creating a lattice-like structure that is most technologically 
useful. This method of analysis can be extended for use in any "nucleationless" model of SAQD 
growth, although here, the specific instance of elastic anisotropy is treated as elastic anisotropy is 
the most easily estimated. Anisotropy of surface energy may also effect SAQD order, and a similar 
analysis would result. 

At later stages when surface fluctuations are large, non-linear dynamics come into play. At 
this stage, there is a natural tendency of SAQDS to either order or ripen.-^ Ripening systems 
will tend to have increased disorder as time progresses, while ordering systems will tend to order 
only slightly due to critical slowing down.- Also, it is important to note that while an ordering 
system might be "ordered" when compared to other non-linear phenomena such as convection 
roles, etc.-, the requirements for technological application are much more stringent. In a numerical 
investigation, it was found that SAQDs have enhanced order due to crystal anisotropy, but soon 
become disordered from ripening.' In any case, an understanding of the order during the initial 
stages of SAQD growth is essential to further investigation of the final SAQD array order. 

As i n'1^1^1^ , a wetting energy is included in the analysis. The wetting potential ensures that 
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growth takes place in the Stranski-Krastanow mode: a 3D unstable growth occurs only after a 
critical layer thickness is achieved, and a residual wetting layer persists. Although somewhat 
controversial, the physical origins and consequences of the wetting potential are discussed in^. 
The analysis presented here is quite general, and one can exclude or neglect the effect of the 
wetting potential by simply setting it to zero. That said, if the wetting potential is real, the present 
analysis shows that it beneficial to SAQD order to grow near the critical layer thickness. 

The remainder of this report is organized as follows. First, stochastic initial conditions are 
discussed. Second, the evolution of a single mode for the isotropic case is discussed. Third, the 
resulting correlation functions and correlation lengths are derived for the isotropic case. Fourth, 
the analysis is repeated for the anisotropic case using elastic anisotropy as an example. Finally, a 
representative numerical example is presented using parameters appropriate for Ge dots grown on 
Si. 

To analyze resulting SAQD order, the mathematical model must include stochastic effects. For 
simplicity, stochastic initial conditions with deterministic time evolution are used. This method of 
analysis yields a correlation function that is used to characterize SAQD order. 

In this model, the film height His a function of lateral position x and time t. The film height is 
treated as an average film height (7Y) with surface fluctuations /i(x, t) , 

n = n + h{^,t). (1) 

In this way, H functions as a control parameter- physically signifying the amount of available ma- 
terial per unit area to form SAQDs, and h(x, t) evolves via surface diffusion giving the the resulting 
surface profile. Order is then analyzed using using the spatial correlation function, (/i(x, t)h(0, t)) 
and the corresponding spectrum function {hk,{t)h]^/{t)*). 

An initially flat surface is in unstable equilibrium, and it is necessary to perturb it to produce 
SAQDs. Therefore, stochastic initial conditions are implemented by letting h{x, 0) in Eq. [T]be 
random white noise. Specifically, h{x, 0) is assumed to be sampled from a normal distribution 
such that 

(/i(x, 0)) = 0, and (/i(x, 0)/i(x', 0)) = A^^'^ (x - x') (2) 

where A is the noise amplitude of dimension [length] ^"^"^Z^, d is the dimension of the surface, 
and 5'^{x) is the rf-dimensional Dirac Delta function. Much of the following analysis uses the 
Fourier transform with the convention, /ik(^) = (27r)~'^ J d'^k exp(— zk ■ x)h{x, t). The mean and 
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two-point correlation functions for /ik(0) are 



(/ik(0)) = 0, and (/ik(0)/ik'(0)*) 



A2 



5'' (k - k') . 



(3) 



The deterministic evolution of a single Fourier component is determined by surface diffusion 
with a diffusion potential yu(x, t) This model is phenomenological in nature, but contains the es- 
sential elements of SAQD formation. Thus, it is an adequate, but not overly complex starting point 
for the investigation of the effects of crystal anisotropy. Furthermore, models of this nature can be 
derived from atomic scale simulations.'" At any instant in time, the growing film is described by 
the curve 7i(x, t). Using Eq.[T]to decompose the film height. 



where fi depends on Ti, and Q is the flux of new material onto the surface. 

The appropriate diffusion potential /i must produce Stranski-Krastanow growth. Thus, it must 
incorporate the elastic strain energy density uj that destabilizes a planar surface, the surface en- 
ergy density 7 that stabilizes planar growth and a wetting potential W{T-C) that ensures substrate 
wetting. The simplest form that gives the appropriate behavior is 



similar to ^^^^^^^ where f2 is the atomic volume, k is the total surface curvature, and is the vertical 
component of the unit surface normal n. The strain energy density u is found using isotropic linear 
plane strain elasticity. In general, u; is a function of x, and it is a non-local functional of the entire 
surface profile H{x). 

Consider first, the the one-dimensional and two-dimensional isotropic cases. FoUowing^ "^ '^, the 
surface diffusion potential, (Eq. Hj) is expanded to first order in the film height fluctuation h. The 
elastic energy u is calculated using linear isotropic elasticity. It is a non-local function of h(x); 
thus, it is useful to work with the Fourier transform. In the isotropic case, the linearized diffusion 
potential (Eq.O depends only on the wave vector magnitude k = ||k|| . Thus, 




dMx,t) = V-(W/x(x,t;7t:)); 



dn/dt = Q 



(4) 




/Uk = f{k,n)hk, with 
f{k,n) = n{-2M{l + u)elk + -fe + w"{n)) . 



(5) 



(6) 
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In the anisotropic case, there will also be a dependence on the wave vector direction 6^. 

The time dependence of the film height has a simple solution if there is no additional flux of 
material {Q = 0, and H is constant). 

hk{t) = /ik(0)e'''=*, with (7) 
ak = -Vef{k,H), (8) 

where ak is the dispersion relation and depends only on the wavevector magnitude k. Modes 
with positive ak are grow unstably, while modes with negative values of decay. 

The important features of are most easily recognized using a characteristic wavenumber is 
kc = 2M(1 + iy)el/-f and a characteristic time tc = {Vn'yk^)~'^r^ 



ak = t-^a^{a- - p). (9) 

where the shorthand a = k/k^ and [3 = W'CH) / {'jk'^) is used. 
The dispersion relation (Eq.[9l) has a peak at ko = a^kc where 

«o = ^ (3 + a/9 - 32/3) . 

Expanding about /cq, 

1 2 

o-fc ~ ctq - {k - ko) where 

f^o = («o - 2/3) , (r2 = k~X\3ao-iP). (10) 

Thus, 

hk{t) ^ /ik(0)e""*~^"2t(fe-feo)^_ (11) 



Now, the statistical correlation functions and correlation lengths that characterize order are 
derived. Using Eqs.|7] and [TT] along with the stochastic initial conditions (Eqs.[2]and[3l) , the mean 
value of hif_{t) is 

(/^k(t)) = (/^k(0))e'^'=* = 0, 

so that the mean surface perturbation is simply for all time and all k. However, the mean-square 
surface perturbations as characterize by the second-order correlation function^^ can be large, 

{h^m^'m = (/ik(0)/ik'(0)*)e('^'=+'^^')* 
5 



using Eq. [3l The real space correlation function can be found by taking the inverse Fourier trans- 
form of Eq. [121 



(/i(x,t)/i(x',t)*) = j d'^^j d'^k' e'^ ''-'^'-''' {hk{t)hk'ity) 

A2 r 

= / d'^ke^^<''~''''^e^'"'' (13) 

(27r)^y 

where integration over k' is simple due to the 5'' (k — k') in Eq.[T2l 
Using EqW 

which is peaked at A; = /cq. This form suggests that the real-space correlation function is periodic 
with a gaussian envelope that has a standard deviation of 

Lcor = V2^ (14) 

Lcor is the correlation length of the film-height profile, and characterizes the degree of order of 
the SAQD array. For example, in one dimension (d = 1), 

{h{x,t)h{0,ty) ... 



^2aot—^Llg^{k—ko)±ikx 



.2 poo 

-y / dk< 



± 

2A2 



an approximation that is valid if kcLcor ^ 1- Thus, Lcor gives the length scale over which dots 
will be ordered, and this scale grows as t^/^. As t ^ oo, 

{h{x,t)hiO,ty) = —-——-e^-^'cos{kox), 

and the entire array should be perfectly ordered. 

The situation in two-dimensions, however, is less friendly. As t — > oo, L^or — > oo, and 
g-iLL.(fc-feo) ~ (2ix/LU^lH{k - ko) 



{h{^,t)h{o,ty) 

A2 



(27r)3/2Leo. 

'^0 ^2aot 



J d^ke^'">'+'^-''5{k-ko) 



A^k 



^^^e-Jo(.o 11x11) 



Thus, the two-dimensional isotropic case has statistical order at large times, but does not yield 
SAQD lattices as does the one-dimensional case (see Fig.[T^ and b). 

Now, consider the effects of crystal anisotropy. For example, let the elastic energy term uj 
have A^-fold symmetry while the other terms are assumed isotropic. Then, the growth rate de- 
pends on the both the wave vector magnitude k (and thus on a = kjk^ and direction 6'^ so that, 

— ^ o"k- Naturally, the anisotropic elastic energy term in the growth rate a-^ depends on the spe- 
cific anisotropic elastic constants, but the general qualitative effect of elastic anisotropy on SAQD 
growth kinetics can be investigated without incorporating a detailed elastic calculation at this time. 
Thus, the present work provides motivation for more detailed calculation. A reasonable way to 
estimate how the elastic energy term varies with direction (6^^) is to assume a low order harmonic 
form with the proper rotational symmetry. The simplest such guess is 

cJk = t-^o? [a (l - esin^ {NB^I2)) - c? - ^) , 

where e parameterizes the importance of the directional-dependence. More precise calculations 
using anisotropic elastic constants of real materials (similar to ref.— ) will be given in future work. 
(Tk has peaks at wave vectors, 

k„ = fco ((cos6'„) 63. + (sin6'„) e^) 

with Qn = 27r(?T, — 1)/A^. Around each peak, k can be decomposed in the direction parallel (/cy) 
and perpendicular (A;^) to k„. Using this decomposition and expanding e^'^'^* about each peak, 

N 

e2-''*^^exp(2(T„t); (15) 

n=l 

lont = 2(rot - ^Lj(A;|j - aok.f - ^Llkl; (16) 
L|| = and = ^^{N^aoet-^k-^)t. (17) 

Note that Ly is the same as Lcor for the isotropic case. Eq.[T5]is valid if k^L^^ ^ 1, and kcL± ^ 1. 
Using, Eq. \T3[ but noting that (Tk now depends on both the magnitude and direction of k, along 
withEqs.[l5l and [IS 



(Mx,t)/i(0,t)*) 



27r 



N/2 
n=l 



(27r)2L||L^ ■■■ 

2 cos {koX\\) 



(18) 
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where x\\ = (cos 9n)x + (sin 9n)y, and x± = (— sin 9n)x + (cos 9n)y. Thus, the same tendency to 
long-range order as for the one-dimensional case is present (see Fig.dh and d). 

As a numerical example, consider Ge grown on Si. Both the isotropic approximation and an 
estimated elastically anisotropic case (with = 4, and e = 0.1) are treated. Neglecting the 
difference in elastic properties of the Si substrate, Eq^. = 1.361 x 10^^ dyne/cm^, z^oe = 0.198, 
eo = -0.0418, n = 2.27 x 10~^\m\ 7 = 1927erg/cm^ and W{n) = 4.75 x 10-V^ erg/cml 
The resulting biaxial modulus is M = 1.697 x 10^^ dyne/cm^, characteristic wavenumber is kc = 
0.369 nm~^ and critical film height is He = 1.132 nm ^ 4 ML. If the film is grown to a thickness 
of = K + 0.25 ML ^ 1.203 nm, and then allowed to evolve, (3 = 0.2078, ao = 0.5664, 
ko = 0.209 nm-\ ctq = 0.01206t-\ ^2 = 0M7k~Hr:\ L,, = 0.746A;(7^(t/t,)i/^ and = 
0.539kQ^{t/tc)^^'^. Note that the unspecified diffusivity V has been absorbed into tc- The film will 
stay in the linear regime as long as the surface fluctuations stay small. For this purpose, let "small" 
mean 1 ML = 2.83 x lO^^cm. Once the fluctuations become "large", individual dots will begin 
to form, and a nonlinear analysis becomes necessary. It is useful to know the correlation lengths 
at this time. 

To find the correlation lengths, one must choose the initial height fluctuation intensity A and 
then calculate the time for fluctuations to become "large". The initial fluctuation intensity is some- 
what arbitrary, but A = 8.02 x 10"^^ cm^ gives an average fluctuation of 1 ML over a patch 1 ML^ 
and seems appropriate. Next Eqs.[l7]and[T8]are used to find t for which the r.m.s. fluctuations 
become large hrms = (1^(0, P)^^^ = ao- There are two solutions, t/tc = 5.53 x 10"'^, and 
t/tc = 471. The first solution is an artifact of the white noise initial conditions and occurs during 
an initial shrinking of the surface height variance; thus, the second solution is taken. At t/tc = 471, 
the correlation lengths are found using Eq. [I7l Lcor = L^^ = 77.5 nm, and L± = 56.0 nm. The 
smaller correlation length gives koL±/7r = 3.73, so a patch of about 4 dots across is expected to 
be reasonably well ordered. 

A numerical simulation of linear size / = 24(27r/A;o) = 722 nm can be easily performed. The 
discrete initial conditions /ik(0) are taken from a normal distribution of zero mean and variance 
(/ik(0)/ik'(0)) discrete ^ ( A^//^)5kk' , whcrc (5kk' is the Kjonccker delta, and cach vcctor compoucut 
of k takes discrete values 2nm./ 1 with m an integer. These components then evolve via Eq. |7J 
The results of an isotropic and anisotropic simulation along with the corresponding real-space 
correlation functions are shown in Fig. \T\ These plots clearly demonstrate the importance of 
anisotropy in producing long range order. 
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Figure 1: Density plots of surface profiles and correlation functions at t/tc = 471. The scale is in nm, and 
the physical parameters are Usted in the text. Plots (a) and (b) are a sample surface profile and correlation 
function respectively for an elastically isotropic material. Plots (c) and (d) are a surface profile and corre- 
lation function for an elastically anisotropic material with 4-fold symmetry. Surface profile plots saturate 
black at a minimum of H = H for clarity. White dots indicate surface peaks. Conelation function plots 
span lb the maximum value / 2. 

Figs, and [It and d appear to agree qualitatively with observations of nucleationless growth 
of Ge^^Sii-a; nanostructures on S i'^'^'^i^^ , although mostly with x < 0.5. Typical observed dot 
arrangements appear to correspond to lower values of P than 0.2078 used for the example as they 
are quasiperiodic but less ordered than Fig.[l]:. Quantitative reporting of correlation lengths would 
assist comparison and possibly enable better tuning of phenomenological models to experiments. 
InP/InGaP nanostructures reported in^^ appear similar. 

From the analysis of the isotropic model, it is clear that long range statistical order (long cor- 
relation lengths) require tight distributions in reciprocal space. This long range statistical order is 
achieved in the large time limit, but this statistical order is insufficient to produce a well-ordered 
array of SAQDs. This lack of usable order is reflected in the real-space two-point correlation func- 
tion of the isotropic model. However, in the one-dimensional and two-dimensional anisotropic 
cases, there is tendency to form a lattice after a long time. In the anisotropic case, there are two 
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correlation lengths that characterize observed SAQD array order. Analytic formulas for these cor- 
relation lengths have been given for a model with elastic simplified anisotropy, but the general 
conclusions and method of analysis should apply to any source of anisotropy. Additionally, from 
a simple form of the wetting potential it is observed that at film heights just above the critical 
threshold Tie, the correlation lengths can grow quickly while the height fluctuations grow slowly; 
thus order is enhanced. Availability of measured SAQD correlation lengths would help to improve 
this analysis and to engineer more ordered quantum dot arrays. 
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